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Abstract 

We analyze the possibility of superluminal neutrino propagation 5t> = (f — c)/c > 
as indicated by OPERA data, in view of previous phenomenological constraints from su- 
pernova SN1987a and gravitational Cerenkov radiation. We argue that the SN1987a data 
rule out 5v ~ (E^/Mn)'^ for N < 2 and exclude, in particular, a Lorentz-invariant in- 
terpretation in terms of a 'conventional' tachyonic neutrino. We present two toy Lorentz- 
violating theoretical models, one a Lifshitz-type fermion model with superluminality de- 
pending quadratically on energy, and the other a Lorentz- violating modification of a mass- 
less Abelian gauge theory with axial-vector couplings to fermions. In the presence of an 
appropriate background field, fermions may propagate superluminally or subluminally, de- 
pending inversely on energy, and on direction. Reconciling OPERA with SN1987a would 
require this background field to depend on location. 



1 Introduction 

Data from the OPERA experiment have recently been interpreted [1] as evidence for su- 
perluminal Ufj^ propagation between CERN and the Gran Sasso laboratory, with 6v = 
{v — c)/c ~ 2.5 X 10~^ for {E,y) ~ 28 GeV ^. Such an extraordinary claim clearly re- 
quires extraordinary standards of proof, notably including confirmation by an independent 
experiment such as MINOS, T2K or NOz/A. Nevertheless, even while the OPERA data 
are undergoing experimental scrutiny, notably of the technical issues of pulse modelling, 
timing and distance measurement on which we are not qualified to comment, it may be 
helpful to present some relevant phenomenological and theoretical observations about the 

^ OPERA used a similar experimental approach to that proposed in [2]. 
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claimed effect. Here we report two sets of considerations concerning: (1) comparison with 

other phcnomcnological constraints on possible superluminal neutrino propagation, and 
(2) instructive theoretical toy models of Lorentz violation that exemplify the price to be 
paid to obtain such an effect. These toy models cast light on possible experimental probes 
of the OPERA effect. 

As we show, reconciling this effect with other bounds on the propagation speeds of 
neutrinos, notably those provided by the supernova SN1987a [3, 4, 2], is a non-trivial 
issue. For example, if 5v were independent of energy, the SN1987a neutrinos would have 
arrived at Earth years before their optical counterparts. This prematurity would have been 
even more pronounced for 'conventional' Lorentz-invariant tachyons, for which 5v would 
increase at lower energies, forcing one to consider Lorentz-violating models. However, 
simple Lorentz-violating power-law modifications of the neutrino propagation speed 5v ~ 
{Ey/Miq)^ arc also severely constrained by SN1987a. Specifically, constraints for = 1,2, 
derived previously in the paper [2] of which one of us (J.E.) was an author, are incompatible 
with the OPERA result for {E^) ~ 28 GeV [1]. Moreover, OPERA reports [1] that there is 
no significant difference between the values of 5v measured for the lower- and higher-energy 
data with {E^) ~ 13 and 43 GeV, respectively, providing no indication that N ^ 0. 

We also discuss the constraints imposed on superluminal neutrino propagation by 
Cerenkov radiation in vacuo. Electromagnetic Cerenkov radiation is suppressed by the 
absence of an electric charge for the neutrino [5]. However, gravitational Cerenkov ra- 
diation [6] is potentially significant for high-energy neutrinos, and an effect of the type 
reported by OPERA could suppress high-energy astrophysical neutrino signals such as 
those associated with the GZK cutoff and with gamma-ray bursters. 

As a complement to these phenomenological remarks, we present two models for Lorentz- 
violating fermion propagation, with different energy dependences for the superluminality. 
One is a simple renormalizable Lifshitz-type fermion model in which the superluminality 
increases quadratically with energy. This model also exhibits dynamical fermion mass gen- 
eration and asymptotic freedom. The other is a modification of an earlier gauge model for 
Lorentz violation with subhimmal neutrino propagation that proposed previously by two 
of us (J. A. and N.E.M.) [7]. The modified model has a background axial U(l) gauge field, 
and may exhibit s?/;?erluminal neutrino propagation with 5v falling as the inverse of the 
energy and depending on direction. Since superluminal propagation with 5v ~ 2.5 x 10~^ 
or greater is incompatible with the SN1987a data, such a model could be compatible with 
the data only if the background field depends on spatial location, and is enhanced in the 
neighbourhood of the Earth compared to its mean value along the line of sight to the 
Larger Magellanic Cloud. 
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2 Phenomenological Constraints on Superluminal Neu- 
trino Propagation 

The primary OPERA result on the mean neutrino propagation speed is 

5v = (2.48 ± 0.28 ± 0.30) X 10-^ (1) 

for z/^ with {EJ) = 28 GeV, where the errors in (1) arc statistical and systematic, respec- 
tively. The OPERA Collaboration also provides the following supplementary information 
on the difference in mean arrival times of samples of higher- and lower-energy neutrinos 
with (E^) = 13 and 43 GeV, At = 14.0 ± 26.2 ns, which corresponds to 

A(5v) = (0.57 ± 1.07) X 10^^ (2) 

for the difference A{6v) between the propagation speeds of these neutrino samples ^. 

Constraints on possible deviations of the speed of z/^ propagation from the velocity of 
light had been placed previously by the MINOS Collaboration [10], which found 

5v = (5.1 ± 2.9) X 10-^ (3) 

for Ufj, with a spectrum peaking at -Ej^ ~ 3 GeV and a tail extending above 100 GeV. The 
MINOS result (3) is not significant in itself, but is also compatible with the OPERA results 
(1, 2), as are earlier neutrino results [8]. 

However, more stringent constraints on models of neutrino propagation are imposed 
by the SN1987a neutrinos [3]. The observed neutrinos emitted by SN1987a had energies 
around three orders of magnitude smaller than the OPERA neutrinos. A significant frac- 
tion of them were undoubtedly z/^, and neutrino oscillation phenomenology severely con- 
strains differences in the propagation speeds of different neutrino flavours, so the OPERA 
results (1, 2) may be confronted directly with the SN1987a data. Since the distance to 
SN1987a was ~ 50 kpc, i.e., ~ 170,000 light-years, an energy-dependent 6v of the mag- 
nitude (1) would have caused the SN1987a neutrinos to have arrived over 4 years before 
their photon counterparts, whereas the maximum tolerable advance is only a few hours, 
corresponding to 5^; ~ 2 x 10~^ [4]. 

The SN1987a data are orders of magnitude more problematic for 'conventional' tachy- 
onic neutrinos. Assuming Lorentz invariance, these would have a dispersion relation 
E"^ = — yU^, where /i^ > 0, and the corresponding deviation of the propagation speed 
from the velocity of light would be 6v ~ ii'^/2E^. Thus, as the energy increases, the speeds 
of such 'conventional' Lorentz-invariant tachyonic neutrinos would decrease towards the ve- 
locity of light. Normalizing the effect to the value (1) of Sv reported by OPERA for their 
relatively high-energy neutrinos would lead to an impossibly large effect for the SN1987a 
neutrinos. Moreover, the magnitude of /i^ would be incompatible with limits on the i/^ 

^Wc note that no other experiment has made such accurate velocity measurements for any other particles 
with Lorentz boosts as large as the OPERA neutrinos. 
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mass from tt and ji decay and, to the extent that oscillation experiments constrain the 
Vf. — mass difference, also the direct constraint on the v^. mass [9]. 

We arc therefore led to consider the possibility of Lorentz violation. Although the 
OPERA Collaboration sees (2) no significant energy dependence of 5v when comparing 
its lower- and higher-energy samples with {E^) ~ 13 and 43 GeV [1], the SN1987a data 
motivate us to look at the implications of an energy dependence 5v ~ {Ej,/M^)^ with 
5v ~ 2.48 X 10-^ for {E^) ~ 28 GeV. Under this hypothesis, the OPERA data would 
correspond to 

Ml ~ 1.1 X 10^ GeV, (4) 
or M2 ~ 5.6 X 10=^ GeV, (5) 

for linear and quadratic energy dependences, respectively. However, stringent constraints 
on Ml and M2 have been imposed previously by observations of the neutrino burst from 
SN1987a [3], which would have been spread out by any energy-dependence of 5v. The 
following constraints on superluminal neutrino propagation were established by a collabo- 
ration including one of the present authors (J.E.) [2] ^: 

Ml ~ 2.5 X 10^° GeV, (6) 
or M2 ~ 4 X 10^ GeV. (7) 

We recall that the supernova neutrino burst is expected to have contained large fractions 
of and v^, so that these constraints apply a priori to the used by OPERA in their 
measurement. 

Comparing (4) with (6), we infer that a linear dependence of 5v is not compatible 
simultaneously with the OPERA and SN1987a data. The situation with a quadratic energy 
dependence, cf, (5) and (7), is not in such stark contradiction with the SN1987a data, but 
the latter would prefer a stronger energy dependence that would be even more difficult to 
reconcile with the lack of any indication of energy dependence within the OPERA data (2). 
As for a possible constant 5v, we recall that the OPERA measurement of 5v ~ 2.5 x 10~^ 
would have led to the SN1987a neutrino signal being observed ~ 4 years before the optical 
signal, whereas the observed advance of < 3 hours (which is compatible with models of 
supernova explosions) would correspond to < 2 x 10~^ [4]. We infer that only an energy 
dependence of 5v ~ E^ with N >2 could reconcile the OPERA and SN1987a data, though 
this is unlikely to be compatible with the lack of a significant energy dependence observed 
within the OPERA energy range (2). 

The possibility of gravitational Ccrenkov radiation has been studied in [6]. The case 
studied there was that of a particle propagating at the speed of light emitting suhluminal 
gravitational radiation, but the same analysis applies to a superluminal particle emitting 
gravitational radiation travelling at the speed of light. It was shown in [6] that a particle 

■^Thc prospective sensitivity of the OPERA experiment to possible superluminal neutrino propagation 
was also estimated in [2], using a similar technique and with with results similar to those now obtained 
(4, 5) by OPERA. 



4 



would lose all its energy within a time 



f - (R) 
tmax - 1)2^3' 

where Mp ~ 1.2 x 10^^ GeV is the Planck mass and n is the refractive index: n = 1 — 5v 
in our case. Setting 5v ~ 2.5 x 10"^ as suggested by OPERA, using (8) we find 

2 X 10^ 

We conclude that applying the OPERA result simple-mindedly would exclude by many 
orders of magnitude the observation of GZK neutrinos [11], which should have ~ 
10^*^ GeV and propagate ~ 10^ light-years ^. Alternatively, neutrinos with ~ 2 x 
10^ GeV, the minimum for which the IceCube experiment has so far published an upper 
limit on the flux [12], could not travel more than ~ 10~^ seconds, ample to explain their 
non-observation, though this surely has a less radical explanation! Conversely, observation 
of neutrinos violating the bound (9) would invalidate the hypothesis of a constant 5v with 
the magnitude suggested by OPERA. 



3 Lorentz- Violating Models with Superluminal Fermion 
Propagation 

In light of the foregoing phenomenological discussion, one might be tempted to lose in- 
terest in theories with superluminal neutrino propagation. However, the effect reported 
by OPERA is so striking and of such potential significance that it is important to study 
whether such an effect is possible, even in principle, and how theoretical possibilities could 
be constrained by future experiments. In this Section, we show how to construct examples 
within the general framework of field theories with higher-order spatial derivatives, and 
discuss some characteristic experimental signatures. 



3.1 Lifshitz-type Field Theory 

Such theories have recently attracted renewed attention because of their improved con- 
vergence properties (for a recent review see [13] and references therein). In this spirit, 
a renormalizable Lifshitz-type theory of gravity has been proposed, which could lead to 

a renormalizable quantum gravity theory at high energies [14]. Such theories arc free of 
ghosts, since the order of the time derivative in the action remains minimal, so that no 
new poles appear in particle propagators. However, in general such theories violate Lorentz 
symmetry at high energies [15]. 

"'As for possible GeV-range neutrinos emitted by gamma-ray bursters with cosmological redshifts, (9) 
shows that they would lose their energy before reaching the Earth, in addition to arriving at very different 
times from their optical counterparts if 6v is given by(l). 



5 



We now exhibit a Lifshitz-type model with superluminal fermion propagation. The 
model is formulated in three space dimensions with anisotropic scaling parameter z = 3. 
In this scenario, the mass dimensions of the coordinates are [t] — —z — —3, [x] — —1, and 
the free fermion action is [16] 

S4ferm = / dtdx (^i^o^ - ^{M^ - A){id ■ t)^' + gi^^f) , (10) 

— * — * 

where A = —did"^ = d ■ d, and we use (+1,-1,-1,-1) as the metric signature. The model is 
renormalizable since [g] = 0, and we have also [M] = 1, [ip] = 3/2. This model exhibits 
asymptotic freedom as well as the dynamical generation of a mass mdyn for the fermion, 
as discussed in more detail in [16], [13]. 

Taking this fermion dynamical mass into account, we obtain the following dispersion 
relation: 

= + MY + 2MY + / , (11) 

and, assuming that M 7^ 0, it is possible to recover approximately Lorentz-invariant kine- 
matics in the infra-red limit, since the rescaling a; = M^a) leads to 

^"-^ln+P' + ^y + ^,, (12) 

where /idyn = "^dj/n/-^^- Using this dispersion relation, one can compute the group velocity 
^ and the phase velocity ^ as power series in [p/MY^p: N > 1. The superluminal 
character of both follows immediately from (12), with the first correction of order p'^/M'^, 
which is quadratically suppressed by the Lorcntz-violating mass scale M. We note that 
the superluminality of this model is an unavoidable consequence of the relative signs of 
the various terms appearing in (10), if one is to avoid tachyonic modes for sufficiently high 
momenta p. 

It is clear from (11) that the superluminality 5v increases quadratically with the fermion 
momentum (or energy) for p{E) < M, and even faster at higher momenta (energies). As 
discussed in the previous Section, such a quadratic dependence is not easy to reconcile 
with the lack of energy dependence in 5v seen in OPERA data [1] , though it comes closer 
to compatibility between the OPERA data and the constraint imposed by SN1987a [2]. 

The fact that the superluminality in this model is quadratic in E^, implies that no effect 
should be seen at the level (1) in the MINOS and T2K experiments, since they have mean 
energies that are almost an order of magnitude lower than the CNGS beam. In particular, 
the indication that Sv from previous MINOS data [10] would not be confirmed in this 
scenario. 



3.2 Lorentz- Violating Gauge Theory 

We now consider more complicated models that lead to forms for 5v with very different 
energy dependences, involving a fermion coupling to either a vector or an axial U(l) gauge 
field. If there is a background field with a suitable constant value in a given reference frame. 
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these models may exhibit superluminal fermions, as well as other dramatic signatures 
highlighted below. These models are concrete realizations of the ideas of [17], where the 
phenomenology of Lorentz violation has been discussed in models where the maximal 
speeds for various particles depend on the species. 

A minimal Lorentz-violating (LV) extension of massless Quantum Electrodynamics 
(QED) was proposed in [18], in which higher-order spatial derivatives were introduced for 
the photon field, and fermions remained minimally coupled to the photon. This theory 
has the features that the light-cone 'seen' by fermions differs from that 'seen' by the 
photon. Specifically, in the theory of [18] (i) the photon always travels at the conventional 
speed of light, (ii) fermions travel subluminally, and (iii) fermion masses may be generated 
dynamically in such a framework, as an alternative to the Higgs mechanism. We will show 
that similar theories with a background vector or axial U(l) field (see also [7]) may lead to 
superluminal fermion propagation, albeit with no mechanism for fermion mass generation. 

The Lagrangians of the models read: 

U,A - -\g,, (l - + ^ (, ^ _ g^^^ ptt) tP-w^tP, (13) 

where G^i, = d^B^ — d^B^ and B^ is a gauge field with either a vector coupling gv or 
an axial coupling qa^ depending whether F = 1 or 75, respectively. The presence of an 
axial 757^* fermion/gauge boson vertex would introduce the possibility of chiral anomalies, 
which could be cancelled by suitable choices of the couplings to the different fermion fields 
■0 = (■^i, • • • , ■0„), represented here by the matrix r with the property: 

tr{r} = . (14) 

In the case of a doublet of fermions, for definiteness, one could use 

although other choices can be made [19], as long as the anomaly-free condition (14) is 
satisfied. It should be noted that no higher-order spatial derivatives are introduced for the 
fermion fields because, in order to respect gauge invariance, such terms would need to be 
of the form 

j^^iiD, ■ 7)> n>2, (16) 

whereby D5 denotes the axial-gauge-field fermion covariant derivative, which would intror- 
duce new, non-renormalizable couplings. The Lorentz-violating modification proposed in 
the Lagrangian (13) does not alter the photon dispersion relation, which remains rela- 
tivistic, but does modify the fermion propagator, as we discuss below. 

It was observed in [20] that models of this type can be obtained by considering the prop- 
agation of photons and charged fermions in a D-particle model of space-time foam [21], 
according to which our world is viewed as a 3-brane propagating in a higher-dimensional 
bulk space that is punctured by point- like DO-brane defects (D-particles). Such models my 
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lead to non-trivial optical properties of the vacuum, because electrically-neutral matter ex- 
citations, such as photons and neutrinos, may acquire non-trivial refractive indices through 
non-trivial interactions with the D-foam. In previous D-foam models [21], these interac- 
tions led to subluminal propagation. In the flat space-time limit, where the low-energy 
Lagrangian is derived, the microscopic reason why fermions do not have higher-derivative 
modifications was that charge conservation forbids interactions of charged fermions with 
the foam. 

The explicit Lorentz violation due to the higher-spatial-derivativc term in the action 
(13) implies that the light-cone 'seen' by the fermions is different from that 'seen' by the 
gauge boson, assumed here to be (almost) massless like the photon, always travels at c, the 
speed of light in vacuo. Specifically, the maximal speed for the fermions is smaller than c, 
as in the vector models of refs. [18, 7]. This may be seen by following the one-loop analysis 
of the fermion wave-function renormalization calculated in [18]. Due to the higher-order 
spatial derivatives in (13), the one-loop quantum corrections to the fermion kinetic terms 
are different for time and space derivatives. A similar computation as the one made in [18] 
yields corrections of the form 

((1 + Zo)ao7° - (1 + Z,)d • 7) , (17) 

where 

Zo = -^Qln(l//x)+ln2-i)+0(/xMn(l//x)), 

Zi = -^Qln(l//.) + |-^ln2)+0(/.Mn(l//.)). (18) 

where av,A = 9va/^'^ ^^'^ = ^/M. We note that the dominant terms in (18), which are 
proportional to ln(l//x), are the same for Zq and Zi. This is to be expected since, in the 
Lorentz-invariant limit: M — )■ 00 and hence // — > for fixed fermion mass, we must have 
Zq — Z\. After redefinition of the bare parameters in the minimal substraction scheme, 
where only the terms proportional to ln(l//i) are absorbed, the fermion dispersion relation 
is 

(1 - ^[2 In 2 - \^ ={\- ^[25/9 - (10/3) In 2])'/ + ^2 _ ^^g^ 
and the fermion phase and group velocities v^., Vg are both subluminal: 

= = 1 - ^ - y ln2) + 0(a^,J < 1 , (20) 

whereas the gauge boson propagates with the standard speed of light in vacuo, c, as 
required by gauge invariance [18]. 

A few remarks are in order at this point. First: in view of (20), the above models con- 
stitute explicit microscopic realizations of the class of Lorentz- violating theories of the type 
considered in [17] , with species-dependent light cones. Secondly, the fact that the constant 
wave function renormalization (18) is found to be less than one, which leads to the sublu- 
minal velocities (20), is a rather general property of quantum field theory, stemming from 
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unitarity [22]. Indeed, in field theories with non-negative-metric states, the wave function 
renormahzation A must satisfy < A < 1, which also implies non- negative anomalous 
dimensions. However, there may be cases, e.g., with derivative interactions [22], in which 
negative anomalous dimensions appear, with the consequence that the wavefunction renor- 
mahzation can be larger than one. It would be interesting to investigate the possibility of 
superluminal propagation in such cases, by analogy with the scenario discussed above. 

We now explore the possibility of superluminal fermion propagation in the context of 
the above theories. To this end, we first consider the possibility of a constant background 
gauge field B'^\ in which case the relevant part of the action (13) reduces to: 

>Cbckgrd ^^{i ^- gv,A ^^^^rr) V - m^^V ■ (21) 

These models fall within the general category of Lorentz- violating extensions of the Stan- 
dard Model [15], as reviewed in the specific case of neutrinos in [23], taking into account 
the available neutrino oscillation data. Depending on the sign of the background field B^'^^ 
in (21), one may have group velocities for the fermions which are superluminal. The quan- 
tum fiuctuations of the axial gauge field would tend to counteract such superluminality, as 
discussed above (20). Nevertheless, for sufficiently weak couplings ay^A and a sufficiently 
strong background field B^^\ the maximal fermion speed may be superluminal. A value 
of 5v ~ 2.5 X 10~^, as reported by OPERA [1], may be arranged in these models with a 
small, perturbative gauge coupling gv,A < 1 and a background field 5*^°^ of appropriate 
magnitude. 

The vector (axial) interaction of (21) has the same (different) signs for left- and right- 
handed fermions, such as neutrinos and their antiparticles which we assume to be 
Majorana fermions. This could lead to a physically important difference between the 
dispersion relations of neutrinos and antineutrinos, and hence apparent CPT violation: 

= \l {P - 9v,aBY + rri^ + Qva^o , 

coj7 = ^J{pT9v,ABy + ^^±9v,ABo ■ (22) 

where the upper (lower) symbols in the combinations ±, =]= refer to the vector (axial) 
case Notice that these dispersion relations are the usual ones for massive particles, 
though with generalised momenta 

Il^ = cu,T9A,vB^, li = PT9A,vB, (23) 

where the upper signs apply to neutrinos, and to antineutrinos with a vector interaction, 
and the lower signs apply to antineutrinos with an axial interaction. 

^In the simple two-flavour axial model (15), the particle of one flavour would exhibit the same dispersion 
relation as the antiparticle of the other flavour. 

^An effect similar to the axial case in (22), but without the flavour structure, could arise purely geo- 
metrically in the propagation of fermions in space-times that break rotational symmetry, such as rotating 
Kerr black holes or axisymmetric Robertson- Walker Universes, as discussed in [24]. Such geometric effects 
stem from the coupling of the spin of the fermions to non-trivial local curvature effects that arise in such 
space-times. 
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— * 

Assuming that the components B, Bo are constants in a local frame of reference, and 

— * 

defining the angle between the three-vectors p and B to be i?, we may write the phase 
velocity following from (22) for high-energy neutrinos with p ^ m as: 

v^f,^^ c^l^^{\B\cos^ - Bo) + ■ ■ ■ , (24) 
P P 

where dots represent higher orders in 1/p. We obtain a similar expression for antineutri- 

— * — # 

nos but with the replacement \B\ — )■ —\B\ and Bq — )■ —Bo in the axial case. However, 
the superluminality associated with (24) does not apply to the group velocity, which is 
subluminal: 

v, = ^ = l- ^MaB' sin^ ^ + m^) + • • • , (25) 

which is the same for neutrinos and antineutrinos, and where the dots represent higher 
orders in 1/p. Note that — 1| is of order 1/p^, unlike the case of the phase velocity, 
where \vph — 1| is of order 1/p 

As a further step, we modify the background space-time in which the neutrino propa- 
gates, exhibiting an extension of this model with superluminal group velocities. We embed 
the model (21) in a modification of Minkowskian space-time with non-diagonal metric 
components that break the rotational symmetry along a specific axis, 

go^ = Vi, i = 1,2,3 (26) 

where \ V\ = V ^ 1 is considered as a small perturbation ^. For constant and homogeneous 
V, the dispersion relations (22) for neutrinos are modified to 

W^W'g^, = , (27) 

where U'^ is given by (23). From this we obtain: 



= -{P- 9v,aB) ■ V + V(p - gv,ABy + m^ + gv,ABo + 0(\/') , 

= -{p- gv,AB) ■ V + ^J{p^gv,ABf + w? ± gv,ABo + OiV^) . (28) 

— * 

Assuming that the components B, Bo are constants in a local frame of reference, considering 

— * — * 

for simplicity the case with \V\ <^ \B\, and defining the angle between the three- vectors 
p and B to be 'd, and that between p and V to be ip, then we observe that Eq. (28) 



^This model has the interesting feature that the (anti)neutrino propagation velocity depends on the 
direction of propagation. This example raises the possibility that, if a constant limiting velocity of light does 
not apply to neutrinos, perhaps the Michelson-Morley experiment should also be revisited for neutrinos? 

®This is motivated by the suggestion that the metric distortion (26) and the axial background case (21) 
may have a common geometric origin, given that they may both be associated with background space- 
time effects, with the vector B pertaining to the coupling of the (anti)neutrino spin to local curvature 
effects [24], as mentioned above. 
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yields the following expressions for the neutrino phase and group velocities for relatively 
high-energies: p = |p| ^ m, \B\: 



Vph — l — Vcosip-\ ■ [\B\cosv — Bq) -\ 

P P 

Vg = 1-V cos^ - ^^'"^ , + OiV^) , (29) 

and similarly for antineutrinos but with the replacement B — > —B in the axial case. Super- 
luminal group velocities of order Sv ~ 2.5 x 10~^, as reported by the OPERA experiment [1] , 
could be obtained for suitable values of the combination —VcosLp > 0. 

The model (21, 26, 29) has several dramatic and testable consequences: 

• The deviation of the neutrino propagation speed from that of light could exhibit 
non-trivial dependence on E^,, due to the combination of terms in (29), that is not a simple 
power law. Thus, compatibihty with the MINOS result [10] is a non-trivial issue, which 
we address below. 

• The neutrino group velocity would depend on the angle of propagation. This means 
that the speed of propagation would, in general, vary sinusoidally during the sidereal day, 
and could even vary between super- and subluminality ^. This modulation would be absent 

— * 

only for V oriented parallel to the Earth's rotational axis. 

• The amount of superluminality would also, in general, depend on the geographical 

— * 

orientation of the neutrino beam. For example, in the hypothetical example in which V is 
oriented parallel to the Earth's rotational axis, the sign of the effect on neutrinos travelling 
northwards (cf, the Fermilab- Soudan neutrino beam) would be opposite to beams travelling 
in a southerly direction (cf, the CNGS neutrino beam), and would be almost null for a 
beam oriented almost east-west (cf, the T2K neutrino beam). Studying the compatibility 
of MINOS data [10] with this model must therefore take into account the ambiguity in the 
orientation of V, as well as the energy dependence of the superluminal effect in this model. 

• It is possible that the orientation and magnitude of V and B vary on an interstellar 
scale, in which case the SN1987a constraint on the neutrino velocity apphes only to an 
average over space and time of the possible superluminality effect, and there is no a priori 
contradiction with the OPERA result. 

• If the neutrino group velocity is superluminal, the corresponding antineutrino group 
velocity in the same direction would also be superluminal in both the axial and vector 
cases. 

Another possibility is that the vector V (26) may be associated with distortions of 
space-time due to the interaction of the neutrino with space-time defects, as in stringy 
D-particle models of space-time foam [21], in which the vector V is associated with the 
average transfer of momentum from the neutrino to space-time defects with which it in- 
teracts during its propagation. In such a case, the metric would be of Finsler type, i.e., 
depending not only on the space-time coordinates but also on momenta. This possibility 

^However, we would not expect any day-night or seasonal dependence, which is consistent with the 
absences of such effects in the OPERA data [1] . 
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is included within our formahsm, but we do not pursue it further here. We note, however, 
that in such models electric charge conservation (which is enforced by gauge invariance) 
prevents charged matter (such as electrons) from interacting non-trivially with the D- 
particle foam [21, 20], so that only neutral excitations (such as photons and neutrinos) 
may be affected by the foam. This may provide a microscopic explanation of the fact 
that for electrons no deviations from special relativity have been observed with a precision 
~ 10-9 [5]. 

4 Summary and Prospects 

The report from OPERA of superluminal neutrino propagation is very surprising, and it 
may well not survive further scrutiny. Moreover, as we have shown in the earlier part of 
this paper, it is subject to constraints from studies of lower-energy neutrinos, specifically 
those emitted by SN1987a [2], and would have implications for higher-energy astrophysical 
neutrinos. In particular, we have argued that the SN1987a data exclude a 'conventional' 
Lorentz-invariant tachyonic neutrino interpretation of the OPERA data. On the other 
hand, as we have shown through the toy models presented in the latter part of this paper, 
it is possible to construct Lorentz- violating theories in which neutrinos travel faster than 
photons, which always travel at c. We have exhibited such models in which the super- 
luminality either increases or decreases with energy. Superluminal neutrinos should not 
be discarded as a phenomenological impossibility, but rather regarded as a scenario to 
be probed and constrained by experiment. In particular, we have shown that the effect 
could depend on the orientation of the neutrino beam. For the moment, the OPERA 
measurement provides a stimulus for investigating such scenarios, but Lorentz-violating 
superluminal fermion propagation should not necessarily be discarded out of hand, even if 
the OPERA result were not to be confirmed. 

Notes added 

A number of papers reacting to the OPERA effect [1] appeared before ours [25]. There is 
some overlap with the phenomenological considerations presented in [2] and here, but the 
models discussed here do not seem to have been discussed yet in this context. 

We also note that, among the extensive literature since our paper was released, it 
has been pointed out [26] that the modified Cerenkov radiation process v ve^e~ is 
potentially an important mechanism for energy loss by superluminal neutrinos. A first 
direct experimental limit on this process and on the distortion of the neutrino energy 
spectrum that it might induce has been reported [27]. We limit ourselves here to noting 
that the rate for this process is very sensitive to the magnitude of 5f , and also to its energy 
dependence. We leave for future work a detailed combined study of the interplay between 
this and other constraints, pending verification of the magnitude of the OPERA result and 
its energy dependence. . 
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In this context, we also note that in any model with general coordinate invariance, such 
as our model (21), where the modified dispersion relations (27) arise as a result of a non- 
trivial metric background, e.g., (26), one may always find coordinate transformations to a 
frame in which the superluminal effects are absent. For the background (26), responsible 
for the superluminal ^-dependent parts of the group velocity (29), such transformations 
are of the Gahlean form i — > — > — VH , which, from the point of view of a passive 
observer, result in a change in the metric Sg^i = —Vi that can cancel the superluminal 
effects in the dispersion relation in that frame. Since the Cerenkov radiation is a physical 
(observer-independent) phenomenon, it cannot depend on the coordinate choice made by 
the observer, whereas the refractive index can, being frame-dependent. Hence, we conclude 
that the arguments of ref. [26] do not apply directly to our second model. We note that 
this argument would imply that, in the transformed frame, the dispersion relations of other 
particles, such as electrons, are affected. However, this is not in contradiction with the 
current bounds for these particles, which are derived in different experimental conditions, 
and specifically in a different reference frame. 
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